CO J=6— 5 in Arp 220: strong effects of dust on high-J CO lines 
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ABSTRACT 

We report new single dish CO J=6-5 line observations for the archetypal Ul- 
tra Luminous Infrared Galaxy (ULIRG) Arp 220 with the James Clerk Maxwell 
Telescope atop Mauna Kea in Hawaii. The J=6-5 line is found to be faint, with 
brightness temperature ratios (6-5)/(l-0), (6-5)/(3-2) of R65/io=0.080±0.017 
and R65/32 =0.082 ±0.019, suggesting very low excitation conditions that cannot 
be reconciled with the warm and very dense molecular gas present in one of the 
most extreme starbursts in the local Universe. We find that an optically thick 
dust continuum, with t{u > 350 GHz) > 1 for the bulk of the warm dust and 
gas in Arp 220, submerges this line to an almost black body curve, reducing its 
flux, and affecting its CO Spectral Line Energy Distribution (SLED) at high fre- 
quencies. This also resolves the C^ line deficiency in this object, first observed 
by ISO: the near absence of that line is a dust optical depth effect, not a dense 
Photodissociation Region (PDR) phenomenon. Finally we briefly comment on 
the possibility of such extreme ISM states in other ULIRGs in the distant Uni- 
verse, and its consequences for the diagnostic utility of high frequency molecular 
and atomic ISM lines in such systems. In the case of Arp 220 we anticipate that 
the now spaceborne Herschel Space Observatory will find faint high-J CO lines 
at u> 690 GHz that would appear as sub-thermally excited with respect to the 
low-J ones as a result of the effects of dust absorption. 
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Subject headings: galaxies: individual (Arp220) — galaxies: ISM — galaxies: 
star burst — ISM: molecules 



1. Introduction 

Since the discovery of its enormous IR luminosity (Lir(8-1000 /xm)~1.6 x 10^^ Lq, Soifer 
et al. 1984; Emerson et al. 1984), representing ~ 99% of its bolometric luminosity, Arp220 
has been a study of extremes when it comes to the ISM conditions in intense starbursts in 
the local Universe. A prominent member in a prominent class of galaxies (e.g. Sanders & 
Ishida 2004), the large IR luminosity and relative proximity (at DL~77Mpc, it is the nearest 
ULIRG), made Arp220 an early target of numerous molecular line observations. Its large 
CO J=l-0 line luminosity revealed a huge molecular gas reservoir (~ 10^" M0) (Young et 
al. 1984; Sanders, Scoville, & Soifer 1991; Solomon et al. 1997), while luminous transitions 
of heavy rotor molecules such as CS J=3-2, and HCN, HCO+ J=l-0 have demonstrated 
that, quite unlike typical quiescent spirals, most of the molecular gas in this spectacular 
merger is very dense, with n > lO^cm"^ (Solomon et al. 1990, 1992). Such studies have 
recently culminated in the most extensive molecular line survey ever conducted for such 
objects (Greve et al. 2009), making Arp220 a galaxy with the best studied molecular gas 
reservoirs in the local Universe. High resolution CO J=l-0, 2-1 intcrfcromctric imaging 
reveals two compact (2r<0.3"; 108 pc) gas concentrations ~1" (360 pc) apart (Scoville, Yun, 
& Bryant 1997; Downes & Solomon 1998; Sakamoto et al. 1999; Eckart, & Downes 2001), 
whose intense starbursts dominate the IR luminosity of the entire system. Finally a major 
advance was recently made with the imaging of the CO J=3-2 emission (and the adjacent 
continua at 860//m) of the warm and dense gas, that also implied substantial dust optical 
depths even at submm wavelengths (Sakamoto et al. 2009). 

We report single dish CO J=6-5 observations of Arp 220 with the James Clerk Maxwell 
Telescope as part of a multi-J CO and HCN line survey of LIRGs (Lir> 10^^ Lq) in the local 
Universe. We demonstrate that an optically thick dust continuum at submm wavelengths 
submerges this CO hue to near-blackbody emission, and strongly affects the emergent CO 
Spectral Line Energy Distribution (SLED) and the C"*" fine-structure line luminosity for this 
extreme starburst. Finally we discuss whether similar ISM conditions are present in dusty 
starbursts in the distant Universe, lowering their observed (high- J)/(low-J) CO ratios well 
below those typical of star-forming molecular gas, and thus limiting their diagnostic utility. 
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2. Observations, line and continuum flux estimates 

We used the upgraded dual-channel W/D (620-710 GHz) receiver at the James Clerk 
Maxwell Telescope (JCMT)0 at 4092 m altitude atop Mauna Kea in Hawaii, operating single 
sideband (SSB), to observe the CO J=6-5 line (z/rest=691.473 GHz) in Arp 220 on March 15th 
2009 under dry conditions (t225ghz~0.04— 0.06). To ensure the flattest baselines possible, but 
also measure the corresponding dust continuum at Aobs=442yum, we used the fastest beam 
switching mode available (continuum mode) with a beam switch frequency of fbmsw=4Hz 
and a throw of 30" (in azimuth). The ACSIS spectrometei0 was used at its widest mode 
of 1.8 GHz (~780kms~^ at 690 GHz) and two separate tunings, yielding an effective band- 
width of 3.235 GHz (~1400kms~^). This was necessary in order to cover the widest known 
CO line in local ULIRGs (FWZl~900-1000kms~^). The typical system temperatures were 
Tsys~(2000— 3000) K, with a median of Tsys~2500K (including atmospheric absorption). 
The beam size at 691 GHz is Ohpbw=8". Good pointing with such narrow beams is crucial 
and was checked every 45mins using both absolute W/D and differential pointing with the 
A3 receiver (230 GHz), yielding o"x~crj/~l-4" (rms) (Figure 1). The aperture efficiency at 
690 GHz is ?7*=0.32, estimated from the Ruze formula for an 80% membrane transmission, 
an rms dish surface accuracy of ~25/im, and a 77*9=0.68 (aperture efficiency for a perfect 
dish with the illumination taper applied at the JCMT), and verified with observations of 
Venus. The flux calibration uncertainty is estimated with repeated observations of compact 
spectral line standards and is ~ 25%. 

Individual spectra were examined, edited for bad channels, and co-added (both W/D 
channels) to yield the flnal spectrum shown in Figure 2 overlaid with CO 3-2 (JCMT), and 
the HCN J=l-0, CS J=3-2 lines obtained with the IRAM 30-m telescope (from Greve et 
al. 2009). The agreement between overall FWZI's and line centers is excellent although the 
CO J=6-5 line becomes weak towards the velocity range where the high density tracer CS 
J=3-2 line (ncr~2x 10^ cm~'^) becomes especially strong (see discusion in section 3). The 
integrated CO J=6-5 line flux is estimated from 

Snne = / S.dV = -?^G(a) / TldV = ^^%^G(a) / T^dV, (1) 
^Av Vs,^^ JAY Va. Jay 

where Fjcmt = 15.62(Jy/K), and 



^The James Clerk Maxwell Telescope is operated by The Joint Astronomy Centre on behalf of the 
Science and Technology Facilities Council of the United Kingdom, the Netherlands Organisation for Scientific 
Research, and the National Research Council of Canada. 

^http: / /www .jach.hawaii.edu/JCMT/spectralJine/Backends/acsis/acsisguide. html 



-4- 







HPBW 



G((7) = l + 81n2 ( — ^ ) , (2) 



with a~o"x~cry the pointing residuals. The aforementioned factor corrects for the flux bias 
resulting when point-like sources are observed with single dish telescopes with given pointing 
rms errors (Condon 2001). This is usually insignificant for modern mm/submm telescopes 
(e.g. at 345 GHz and the same pointing rms error: G~1.064), yet it can become signifi- 
cant at very high frequencies, even for the excellent tracking and pointing of the enclosed 
JCMT. For our observations G=1.17, and we use it to scale our measured CO J=6-5 line 
and 442/im dust continuum fluxes accordingly. We obtain Sco(6— 5) = (1170±341) Jykms"^, 
and S442^j^=(3.71±0.96) Jy (using the line-free part of the spectrum and the scaling fac- 
tor in Equation 1), or equivalently S434^jjj=(3.85±0.99) Jy (scaled by Sj^ocz/^). Regard- 
ing the continuum it is worth mentioning that Arp220 is a submm calibration source for 
SCUBA at the JCMlfl, where multiple images show a compact source at 450/im (< 5") with 
S450/xm=(2.77±0.06) Jy. The higher value of S45o^m=(6.3±0.8) Jy reported by Dunne & Eales 
(2001) is thus likely affected by calibration error (Dunne 2009, private communication). Fi- 
nally the CO J=6-5 flux is in good agreement with that reported recently by Matsushita et 
al. 2009 of Sco(6-5) = (1250±250)Jykms-^ using the SMA0. For the purposes of this work we 
adopt: Sco(6-5) = (1210±240) Jykms-i (average: JCMT and SMA), Sg4^^=(3.0±0.06) Jy 
(average: W/D, SMA, SCUBA cahbrator database). 



3. The state of molecular gas and dust in Arp 220 

A large multi-J CO, HCN, HCO+ and CS line survey by Greve et al. (2009) makes 
Arp 220 the ULIRG with the best-studied molecular gas properties, and allows us to place 
the measured CO J=6-5 line in the best possible perspective. From that study it be- 
comes apparent that densities for the bulk of the molecular gas in Arp 220 are very high: 
n(H2)>10^^^ cm^^, and thus the CO J=6-5 transition is expected to be fully thermalized 
(ncrit~ 5.8x 10^ cm~^). On the other hand global dust continuum SEDs (e.g. Lisenfeld, 
Isaak, & Hills 2000), and high resolution mm/submm imaging (e.g. Downes & Eckart 2007; 
Sakamoto et al. 2008) yield typically Tdust~ (65-120) K. For the high gas densities in this sys- 



http://www.jach. hawaii.edu/JCMT/continuum/calibration/sens/potentialcalibrators. html 

fringe pattern present on the short baselines has likely corrupted the SMA channel map spectrum 
(Fig. 7 in Matsushita et al. 2009), and thus the SMA spectrum is not reproduced here. The reported SMA 
line flux has been corrected for that effect (Matsushita 2009, private communication) 
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tern: Tkin~Tdust, and the CO line emission should be peaking at Jmax~[2kBTkin/Eo] ~5-7 
(Eo/kB~5.5 K). Clearly the J=6-5 transition should not only be luminous but quite possi- 
bly the most luminous CO line in this system, with typical brightness temperature ratios of 
R'65/j+i,j~0.80-0.95 for J-|-l=l-5, even for the lowest possible Tk=65 K and n(H2)=10^ cm~^. 



Our measured CO (6-5) / (1-0) and (6-5) / (3-2) brightness temperature ratios for Arp 220 
are: R65/io=0.080±0.017 and R65/32 =0.082 ±0.019 (CO J=l-0, 3-2 fluxes from Greve et al 
2009), suggesting impossibly low excitation conditions, very much incompatible with the 
average properties of its molecular gas reservoir. High dust optical depths at short submm 
wavelengths provide the most plausible explanation for the faint CO J=6-5 line, immers- 
ing it in a strong nearly blackbody (and thus almost featureless) dust continuum. Such 
large optical depths of dust continuum at far-IR wavelengths were first proposed for Arp 220 
by Lisenfeld et al. (2000) where Tioo~5-12 has been deduced from its far-IR/submm dust 
SED, with Gonzalez- Alfonso et al. (2004) arriving at similar conclusions using a full ISO 
LWS spectrum. This has culminated with the recent work of Sakamoto et al. 2008 where 
large optical depths are reported even at submm wavelengths (rgeo^l)- Their effect on the 
emergent CO line ratios can be easily shown using a simple model of an isothermal mixture 
of molecular gas and dust in LTE, where the continuum-subtracted source line brightness 
temperature is (e.g. Rohlfs & Wilson 1996): 



where J(i/, T)=hz//kB [exp(hz//kBT) — 1]~ , and ra and Tco are the dust continuum and CO 
line optical depths. Here is important to point out that the equation above assumes a 
form identical to that for line absorption by a foreground dust screen only for an isothermal 
mixture of line-emitting gas and continuum-emitting dust (i.e. identical source functions). 

In such a mixture the high dust optical depths do not actually reduce the emergent line 
strength (as they would for a dust absorption "screen") rather than make both continuum 
and spectral line emission rise up to a common blackbody SED, with Equation 3 simply 
expressing the diminishing line-continuum contrast. For a dust cmissivity law Ta[u)(xu^ 
(/3=l-2), the brightness temperature ratio CO(6-5)/(J-|-l-J) transition would then be 



3.1. High dust optical depths at submm wavelengths 



AT, 



b,co 



e-^MJ(i^co,T)-J(l/eo,TcMB)] (1 



e , 



(3) 
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where R, 



(int) 

'65/J+l,J 



is the "intrinsic" line ratio corresponding to Td—0. For CO J=3-2 (also 



tracing the warm star forming molecular gas as the J=6-5 line) and its corresponding con- 
tinuum emission Sakamoto et al. (2008) has deduced T(^(z>'32)~1- For n(H2)=10^ cm^^ and 
a typical Tkin=(65-90) K our Large Velocity Gradient (LVG) code yields Rg™*32~0-85-0.90. 

Inserting these values into Equation 4 yields Rggy32~0. 042-0. 32 (for /3=l-2), comfortably 
encompassing the observed value, and demonstrating how substantial submm dust optical 
depths can suppress (high-J)/(low-J) CO hne ratios. The larger suppression of CO J=6-5 
in the eastern nucleus of Arp220, contributing mostly to the "horn" at ~5600kms~^ of the 
CO line profiles, is then a likely result of the larger gas densities (and thus correspondingly 
larger dust optical depths). This is clearly indicated by the rising CS J=3-2 hne emission 
(Figure 2) and the larger HCN(4-3)/CO(3-2) ratio (Figure 3) towards ~(5500-5600) kms"^ 
than towards ~ (5300-5400) kms~^, the latter range associated with the western nucleus. 

Thus the CO J=6-5 line (and partially even J=3-2) is affected by a dust continuum 
rising to an almost black body curve from IR out to short submm wavelengths. This 
is further indicated by the observed ratio S434juin/S860/um~5.45±0.82 (for the two nuclei: 
S860Mm=(0.55±0.082) Jy, from Sakamoto et al. 2008), which is close to the black body value 
of (860/434)2=3.93. Indeed from 



and S434^ni/S860Mm=5.45 we find T86o~1-25 (for /3=2), while for S434^in/S860Mm=6.27 {+a of 
the measured value) the minimum optical depth is rggo ~0.95. Dust optical depths large 
enough to bring all this about are certainly present in Arp220 (Sakamoto et al. 2008), while 
the globally supressed J=6-5 line further demonstrates that much of the large molecular gas 
mass in this system is "cloaked" by these high dust optical depths rather than only small 
sub- regions around e.g. an AGN (Downes & Eckart 2007). If this is typical for ULIRGs, it 
can modify the emergent CO SLEDs, and affect the diagnostic value of molecular lines at 
high frequencies for such extreme systems, a possibility further discussed in Section 4. 

Finally we note that the average state of the molecular gas in Arp 220, as indicated by 
the HCN, CS, and HCO"*" lines, is so highly excited that even CO J=6-5 line fluxes up to 
~3 times higher would leave R65/32 well below what is expected of such gas and is typically 
observed in starbursts (see Figure 5 in section 4.2), still suggesting substantial dust optical 
depths at short submm wavelengths surpressing high-J CO line emission. 
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3.2. Arp220: the C+ line luminosity deficit resolved 

For a Tgso ^ 1 "cloaking" the bulk of the gas and dust in Arp220, and conservatively 
assuming /3=1, yields risg ^ 5.4 at the emission wavelength of the C+ fine structure line. 
This is more than enough to almost completely suppress this strong ISM cooling line into a 
featureless black body dust continuum, favoring large dust optical depths as the cause for 
its weakness among the other explanations proposed (e.g. Malhotra et al. 1997). Indeed 
the faint CO J=6-5 fine, besides corroborating the high dust optical depths at IR/submm 
wavelengths in Arp220, also makes them the most likely cause of its C"*" line luminosity 
"deficit" (found in ULIRGs, with Arp220 having the largest; Luhman et al. 1998, 2003). It 
does so by negating a prominent alternative explanation for the line surpression, namely 
very dense PDRs (where precipitous C"*" recombination would remove it from the ISM). This 
is because CO J=6-5 is a different spectral line (governed by different physics and chemistry 
than C"*") and one that ought to be luminous if that alternative explanation held. Indeed, 
for dense PDRs immersed in strong far-UV radiation fields high-J CO lines are expected to 
be very luminous and an alternative cooling "channel" to that of the suppressed C"^ line, 
balancing the tremendeous ISM heating expected in ULIRGs (Papadopoulos, Isaak, & van 
der Werf 2007). Finally it must be noted that if high dust optical depths at short submm 
wavelengths are responsible for surpressing both the high-J CO and the C"*" lines in ULIRGs, 
it follows that the starburst systems with dust-surpressed (high- J)/(low-J) CO line ratios 
will also be those with small L(C+)/Lir. 



4. High dust optical depths in ULIRGs: from the IR to the submm 

The discovery of dust continuum optical depths that in ULIRGs can remain substantial 
even out to short submm wavelengths is the culmination of a series of studies successively 
"pushing" the Tx > 1 limit longwards in wavelength for these remarkable systems (Condon 
et al 1991; Solomon et al. 1997; Lisenfeld et al. 2000; Sakamoto et al. 2008). Condon et 
al. were the first to point out that high dust extinction at far-IR wavelengths is needed to 
explain radio continuum sizes of ULIRGs that are smaller than their minimum IR black- 
body emission sizes. They argued that re-radiation of the IR light from a compact starburst 
by dust "layers" further out that remain optically thick at far-IR wavelengths can yield 
larger effective IR-continuum source sizes than the true starburst sizes (revealed by radio 
continuum). For the compact gas disk configurations in ULIRGs (expected in mergers and 
their dissipative gas motions) such large dust optical depths can make both IR "colors" and 
ISM line ratios at short submm wavelengths viewing-angle dependent and thus much reduce 
their AG N-vs- starburst diagnostic power. This is indeed the case for IR colors where "IR- 



-8- 



cool" ULIRGs may not be lacking an AGN but instead host AGN that are heavily obscured 
along the particular line of sight observed (Condon et al. 1991), while in nearly face-on 
ULIRG/QSO systems (e.g. Mrk231) such AGN remain visible with their warm IR colors 
discernible through smaller dust optical depths. It now seems that molecular line diagnostics 
at short submm wavelengths may suffer from similar effects. 

4.1. The diagnostic power of fEir-IR/submm ISM lines revisited 

A general effect of high dust optical depths at short submm wavelengths on the CO 
SLEDs would be to make them appear "cooler" than those typical for star-forming gas. 
Thus dust-suppressed CO (high- J) / (low- J) ratios can be hard to discern from those typical 
of genuinely low-excitation gas. High dust optical depths at short submm wavelengths in 
ULIRGs are particularly troubling for ISM lines whose physics makes them diagnostic of 
deeply "buried" AGN versus starbursts as the energy source of their prodigious IR luminosi- 
ties. This is especially true for line frequencies oiu > 690 GHz whose ratios, when unaffected 
by dust absorption, offer the best discrimination between the various ISM excitation mech- 
anisms, such as X-rays from a deeply buried AGN (Meijerink & Spaans 2005), or far-UV 
photons from their starbursts (Meijerink, Spaans, & Israel 2006). In such cases ratios of 
lines adjacent in frequency (e.g. HCO^/HCN for the same rotational level) will have to be 
used as AGN/starburst discriminators. Finally, high dust optical depths at short submm 
wavelengths can affect molecular line observations of ULIRGs at high frequencies with the 
upcoming Herschel Space Observatory. For Arp 220 in particular faint and apparently sub- 
thermally excited CO lines are to be expected at i/>690GHz (Figure 4), and this will also 
affect high- J transitions of any molecules tracing the star-forming gas phase (e.g. HCN). 

4.2. Evidence for high dust optical depths at submm wavelengths in distant 

steirbursts 

Low (high-J) / (low- J) CO line ratios quite atypical of star-forming molecular gas have 
been measured in some of the most IR-luminous (Lir~10^^ Lq) starbursts in the distant Uni- 
verse (Tacconi et al. 2006), a Ly- break galaxy whose CO(7-6)/(3-2) brightness temperature 
ratio is R76/32=0.030 (Baker et al. 2004), and in a submm-bright ULIRG at z—lAA where 
R54/2i~0.16 (Papadopoulos & Ivison 2002). However the wealth of molecular line data for 
Arp 220 that allowed its low CO (6-5) / (3-2) ratio to be attributed to large dust optical 
depths at short submm wavelengths is not available yet for high rcdshift systems. Thus, as 
discussed in the previous section, "cool" CO line ratios can be due to either low-excitation 
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gas or to high dust optical depths. 

Dominant amounts of low-excitation and SF-idle gas extending well beyond star-forming 
galactic nuclei are indeed well-known features of star-forming spirals in the local Universe 
and even of archetypal starbursts such as M82 (Weiss, Walter, & Scoville 2004). However 
recent high resolution CO imaging of submm-bright galaxies (SMGs) at high redshifts finds 
no evidence for such low-excitation extended molecular gas reservoirs (Tacconi et al. 2006; 
lono et al. 2009), yet it also measures CO (6-5)/(3-2), (7-6)/(3-2) line ratios of < 0.35 in 
many such systems. In one extreme such case, the most compact SMC (SMM J044307-I-0210 
at z=2.5) has the lowest such line ratio of R76/32~0.13 (Tacconi et al. 2006). Indicatively, for 
nearby star-forming galactic nuclei such CO line ratios are ~0.8-0.9 (Mao et al. 2000; Bayet 
et al. 2006), and values of <0.30 are typical of low-excitation SF-quiescent gas. In Figure 5 
we show the R65/32 ratios for star-forming and quiescent environments in the local Universe 
along with those available for systems at high redshifts. From this it becomes obvious that: 
a) nearby starburst environments behave as expected, tracing the entire range of R65/32 values 
typical for star-forming molecular gas, and b) several high luminosity SMGs (Lir>1O^^L0) 
fall well below this rangejfl and well into that typical for low-excitation SF-idle environments. 
In some SMGs R65/32 falls even below the minimum value of ~0.3, expected for the dense 
and cold starless cores with minimum temperatures of Tkin~10K (set by Galactic cosmic 
rays), and densities of n(Il2)~10^ cm~^. Indicatively R65/32 for Arp220 would remain below 
that value even for a CO J=6-5 line flux that is ~3 times higher than the one observed. 

In order to better quantify this we conduct Large Velocity Gradient (LVG) radiative 
transfer modeling of the CO line ratios reported for SMGs by Tacconi et al. (2008). For 
R65/32 < 0.3 and R76/32 < 0.15 we find no solutions within Tkin~(40-90) K, n(H2) >10^ cm'^ 
and Kvir~l (self-gravitating gas phase), the parameter space typical for star- forming gas. 
The Kvir parameter is 



^ (dV/dr)LVG . .4ICO/H2 
(dV/dr)„^,^, ~ • y^A,o 




(a~l-2.5, depending on the cloud density profile) and quantifies the kinematics of the gas 
responsible for the CO emission, with Kvir~l (within factors of ~2-3) indicating a virialized 
gas phase and Kvir^l an unbound one (e.g. Greve et al. 2009). Only for R76/32~0.3 are there 
solutions Tkm=[40K, (65-80) K], n(H2)=10^ cm~^, i.e. typical for star-forming gas, except 



^If a "neighboring" CO line ratio (e.g. R76/32) rather than R65/32 is avalaible for a high-z system, we 
obtain the average value and range of Rgs /32 by using the available one to find compatible LVG solutions 
over Tk=(30-110)K, (n(H2)) = (102-10^) cm'^, assuming also Kvir-'l. 
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that Kvir=[4, 11] indicates unbound gas motions. The low CO (high- J)/(low- J) hne ratios 
found in some distant starbursts are thus not compatible with the typical conditions of warm 
and dense star-forming gas. Hence, unless one is willing to postulate dominant amounts of 
cold and/or diffuse gas in some of the most spectacular starbursts in the Universe, high dust 
optical depths at rest-frame submm wavelengths are the only other way of surpressing such 
(high- J)/(low- J) CO line ratios. Finally, it is unlikely Arp220 will be the only ULIRG, near 
or far, where this happens making its CO SLED appear "cool" and atypical of its vigorously 
star-forming dense gas. 

In summary high dust optical depths at short submm wavelengths in extreme starbursts 
could then seriously impact 

• the detectability of CO J-l-1— >J, J-|-l>6 lines of starbursts in the local or distant 
Universe with Herschel and ALMA respectively 

• the diagnostic power of high-frequency molecular and atomic line ratios, (e.g. deduced 
densities and temperatures, AGN-induced XDRs versus star burst-related PDRs) 

• the thermal balance in the molecular ISM of ULlRGs (with C^ and high-J CO line 
cooling diminished, dust continuum maybe the dominant cooling "channel" ) 

We will further explore these issues in a future paper that will include the entire suite of 
CO hnes observed for the LIRGs in our sample and their corresponding molecular SLEDs 
for heavily dust obscured environments. 

5. Conclusions 

We report measurements of the high-excitation CO J=6-5 line and its adjacent dust 
continuum at 442/im of the prototypical Ultra Luminous Infrared Galaxy Arp 220 with the 
James Clerk Maxwell Telescope in Hawaii. Our findings can be summarised as follows: 

1. The CO J=6-5 line is very faint, with CO (6-5)/(l-0), (6-5)/(3-2) brightness tem- 
perature ratios of R65/io~R65/32~0.08, much lower than those expected from the warm 
and dense star- forming molecular gas that is present in this extreme star burst. 

2. High dust optical depths with t(i/ > 350GHz)>l are responsible for the faint CO 
J=6-5 line in Arp 220 and will strongly modify its emergent CO SLED. We anticipate 
that for this ULIRG the now spaceborne Herschel Space Observatory will observe faint 
and apparently sub-thermally excited CO lines beyond ~690 GHz. 



- 11 - 



3. The C"^ and the CO J=6-5 hne luminosity "deficits" in this system have the same 
cause: high far-IR/submm dust optical depths. 

4. The low CO (high-J) / (low-J) line ratios found also in several starbursts at high redshifts 
could also be due to high dust optical depths at short sub-mm wavelengths, though 
in individual objects this can be hard to distinguish from large reservoirs of genuinely 
low-excitation molecular gas. 

5. High dust optical depths at far-IR/submm wavelengths can affect the diagnostic power 
of molecular and atomic lines in the spectral regime where it is expected to be the 
greatest, making their intensity ratios depend on differential absorption and viewing 
angle rather than underlying gas excitation conditions. 
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NOTE-ADDED-IN-PROOF: Recent results from Herschel Space Observatory suggest 
a CO J=6-5 line-integrated flux 2-2.5 times higher than that measured by the JCMT and 
the SMA. Pending solution of various calibration issues with the Herschel FTS instrument, 
we note that such values still leave our main conclusions intact (section 3.1) 
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Fig. 1. — The distributions of pointing offsets obtained during all W/D observing periods. 
Values with |cr|<2" are those that apply to source observations, with any larger ones typically 
obtained after large changes in azimuth and elevation when the telescope changed sky sectors. 
The pointing offsets during Arp 220 observations are: aaz~crez~1.4", and Or={o1^+al-^)^^'^r^2" . 
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Fig. 2.— Arp220: a=15^ 34°' 57.24^ 5=+23° 30' 11.2" (J2000). Top: CO J=6-5 (thick 
line), overlaid to CO J=3-2 (thin line), with resolutions of AVch=27kms~^ (J=6-5) and 
AVch=23kms^^ (J=3-2), and centered at cz=5450kms~^ (LSR). The thermal rms error 
across the line-free part of the spectrum is 5T^~12mK for both lines. Bottom: the CO 
J=6-5 line overlaid with HCN J=l-0, and CS J=3-2 lines, obtained with the IRAM 30-m 
telescope (adapted from Greve et al. 2009). 
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Fig. 3. — The HCN J=4-3 (thin line) overlaid to CO J=3-2 (thick line) spectrum, with 
resolutions of AVch=22-23 km s"^ (HCN J=4-3, CO J=3-2). The CO and HCN emission of 
the denser eastern nucleus corresponds to the feature centered at ~5500-5600 km s"^ where 
the HCN/CO line ratio is the largest (Greve et al. 2009). The thermal rms errors across the 
line-free part of the spectra are: (5TX~2mK (HCN J=4-3) and 5TX~12mK (CO 3-2). 
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Fig. 4. — The CO SLED of Arp220 and the effects of dust (see Section 4.1), where Lune 
is given in solar luminosities. The red bars denote the 5sigma-lhr detection limits for the 
SPIRE-FTS on Herschel, evaluated using the observatory integration time calculator and 
assuming that all the line flux is contained in 1 (CO J=4-3 up to J=7-6), 2 (CO J=8-7 up 
to J=10-9) and 3 resolution elements (CO J=ll-10). It is clear that CO lines above J=6-5, 
7-6 can be faint and appear as sub-thermally excited, while very high-J transitions such as 
J=10-9, 11-10 will be undetectable. 
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Fig. 5. — The CO (6-5)/(3-2) ratio versus LiR(40-400/i.m) for: a) quiescent systems (dia- 
monds), nearby starburst galaxies or starburst nuclei (triangles), and c) SMGs (stars). The 
ULIRGs Arp220, Mrk231, the star- forming region in Orion A, and the SF-quiescent Galactic 
Center are also shown. The shaded areas mark the R65/32 range for typical star-forming gas 
in LIRGs (hatched), and the extreme star-forming environments of ULIRGs (cross-hatched). 
The first was estimated from LVG models for: Tkin=(30-100) K, n(H2) = (10^-10^) cm'^ and 
Kvir=l (fully encompassing the typical conditions of star- forming gas), while the second 
narrower range is obtained from LVG solutions constrained by the HCN(3-2)/(l-0) aver- 
age brightness temperature ratio of (R32/io(HCN))=0.55 of nearby ULIRGs (from Krips et 
al. 2008; Gracia-Carpio et a. 2008). The dotted line marks the lowest possible R65/32 for 
dense and cold CMC cores (n(H2)>10^cm^'\ Tkin=10K) where both transitions would be 
fully thermalized and optically thick, and which is also the limit below which only gas with 
n(H2)<10^cm~^ can be found. CO data: Galactic Center: Fixsen et al. 1999; Orion A: Mar- 
rone et al. 2004; Mrk231: Papadopoulos et al. 2007; Arp220: this work; SMGs: Tacconi 
et al. 2006, Papadopoulos & Ivison 2002; nearby starbursts: Wild et al. 1992, Mao et al. 
2000, Bradford et al. 2003, Bayet et al. 2006. 



